Dark matter equation of state from rotational curves of galaxies 
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We obtain a dark matter equation of state by using the rotation curves of galaxies only. Dark 
matter is modeled as a fluid component with pressure and we show that, within the Newtonian 
limit, the resulting equation of state is compatible with existing cosmological bounds. An upper 
bound on the central pressure of the dark matter is obtained, within the full general relativistic 
treatment, by demanding that the dark matter fulfills the positive energy condition. 
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I. INTRODUCTION 

Astrophysical observations of the movement of stars 
in galaxies, of the cosmic microwave background (CMB) 
anisotropies, the weak gravitational lensing of distant 
galaxies, among others, are all explained by the presence 
of dark matter, which has not shown any interaction with 
the baryonic matter except the gravitational one [l[ . On 
the other hand, direct detection of dark matter demands 
it to have some interaction with the rest of particles of 
the standard model 0. These two research lines seem 
to have contradictory points of view on the properties 
of the dark matter. In order to attempt to conciliate 
these seemingly opposite strategies, one can research on 
the possibility that the astrophysical description of dark 
matter could relax the pressure-less hypothesis, and al- 
low some interaction with itself and with the baryonic 
matter. 

For instance, at the cosmological level, Miiller in 
3 proposed a barotropic relation for the dark fluid, 
p = oj-qm Pi instead of considering the usual dust case 
(wdm =0). It was shown that in the case where the adi- 
abatic sound speed coincides with the sound speed, the 
combined analysis of data from anisotropies of the CMB, 
supcrnovac la and matter power spectra strongly con- 
strain udm to lay within the interval [—1.5, 1.13] x 10~ 6 
(see also 0, 

Since the strongest evidence for dark matter at galac- 
tic scales comes from the rotational curves, it is natural 
to look for a functional relation between pressure and 
energy density by modeling the dark matter halo as a 
perfect fluid. In the present work we show that indeed 
it is possible to derive (not to propose) such functional 
relation solely with the use of the rotational curve pro- 
files. We performed a general relativistic treatment of 
the problem, as well as its Newtonian counterpart. The 
relativistic result gives a maximum value of the pressure 
at the center of the galaxy, assuming that dark matter 
should fulfill the positive energy condition. The general 
relativistic result contains, as expected, the Newtonian 
description as a particular case, and, in this limit, the 
derived equation of state is consistent with the bound 



obtained by Miiller In fact, when the value of the 
pressure at the center is very small, there are no remark- 
able differences between both treatments, and the New- 
tonian model can be used with confidence. We found 
an analytical expression for the dark matter equation of 
state in the Newtonian limit as a function of the asymp- 
totic circular velocity of the halo. 

We have assumed spherical symmetry and a static and 
isotropic description for the dark matter fluid. These 
considerations reduce the number of unknowns, so that 
only one observation (the rotational curve profile) is 
enough to determine the equation of state. If the 
isotropic assumption is relaxed, then it is needed another 
observation to determine the equation of state. Such ob- 
servation can be the gravitational lensing, see for instance 



The manuscript is organized as follows: In section |TT] 
we present the general relativity equations for a perfect 
fluid and discuss how one of the metric coefficients is 
completely determined by the rotational curve profile; we 
also present the Newtonian treatment. In section Hill we 
compute the equation of state by means of the rotational 
curves of LSB galaxies fitted by the Salucci profile |{|. 
We finally discuss some of the implications of the derived 
equation of state and conclude in section IIVI 



II. PERFECT FLUID HALO 

We consider a static and spherically symmetric space- 
time in General Relativity, described by the line element: 
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sin 2 9dip 2 . The gravitational potential 



where dfl = dO 2 
$(r) and the mass function m(r) are functions of the 
radial coordinate only. 

One of the geometric potentials can be determined by 
the observations. Indeed, it can be obtained an expres- 
sion relating the tangential velocities of test particles, Vt, 
following stable circular orbits in this spacetime, with 
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the gravitational potential $, by means of solving the 
geodesic equations for test particles in circular stable or- 
bits (see [10l - [l2| for details on the derivation). Such ex- 
pression is 



with 
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where /3 2 = ^. 

It is remarkable that the same expression is obtained 
by the Newtonian description, by simply equating the 
gravitational force with the centrifugal one for particles 
in circular orbit. 

It is convenient to define dimcnsionlcss variables for 
the mass and the radial distance, namely m = Mo n 
and r = R x, where M ,R are the characteristic scales 
for mass and distance of the system under study, and 
n is a dimensionless function related to the mass and x 
a dimensionless quantity related to the distance. With 
these scales, we also define a characteristic density, po = 
Mo/(4/3)7ri?o 3 , and a characteristic pressure: po = poc 2 , 
so that p — pop, and p = po p, with p, p dimensionless 
functions related to the density and the pressure respec- 
tively. 

With these definitions, we rewrite the Einstein's equa- 
tions and the conservation equation as follows: 
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where we have already used the relation ([2]) and prime 
stands for derivative with respect to x. Also, it was de- 
fined q = Mo G/c 2 Rq which is a dimcnsionlcss quantity. 

Thus, to solve the system of a self- gravitating isotropic 
fluid, we are left with three unknowns: one geometric, 
the mass function, and two related to the dark fluid, the 
pressure and the energy density. 

From observations, it is known that the terminal ve- 
locity is in general of the order of 100 Km/s, so that, 
using the Newtonian relation for the mass of the halo: 

m(r) = rv t 2 /G, and taking G = 4.299 10" 6 Kpc £^ B) ' , 
and the characteristic distance in Kpc, we see that the 
characteristic mass of these systems is 10 10 M©. Also, we 
find that the parameter q is of the order of 5 x 10~ 7 . 
We will use this value for q in our computations. With 
these values, the characteristic values of the density and 
pressure of the system are: po = 1.6210~ 22 |^§ and 
Po = 1.46 lO- 1 ^. 

1 v cm s z 

After some manipulation of equations (0[5|), we obtain 
an equation involving only the pressure, p: 



S(x) = 
T(x) = 
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The general solution for the pressure, expressed in 
terms of the observed rotational velocity profile, /3(x), 
has the form 
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where the value of the integration constant, C, is set by 
the appropriate boundary conditions. 

The density and the mass of the dark perfect fluid can 
be directly computed from Eqs. ([3][5|), once the pressure 
is determined: 
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Once we have presented the general rclativistic treat- 
ment for the dark perfect fluid, it is interesting to have 
a look on the Newtonian description in order to asses 
its accuracy. Considering an static spherical dark matter 
distribution we have the Poisson's equation and the defi- 
nition of the mass. To these equations, we add the Euler's 
equation in this case, for the co-movil observer. Thus, we 
obtain the following set of Newtonian equations: 
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p' + S(x)p = T(x), 



(6) 



where the considerations for the defining dimensionless 
quantities are as those given above; we used the relation 
valid in the Newtonian limit, between the mass and the 
gravitational potential, $' = Gm/r 2 , and we express this 
last one in terms of the observed velocity profile ratio, 
p{x). 



III. DARK MATTER EQUATION OF STATE 
FROM ROTATIONAL CURVES 

The dark matter equation of state is directly computed 
from Eqs. ([5]) and © once the velocity profile /3 2 (x) is 
given. We will consider an observational velocity profile 
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proposed by @, which was obtained by adjusting 1023 
galaxies. It is given as: 



/3 2 = Po 



(14) 



be fitted to observational data. 

It is worth to mention that one of the integrals for 
the pressure in Eq. ([8]). is solved analytically when the 
rotational curve profile Eq. (fT4")l is used. Indeed: 



where pb,a are constants, po stands for the ratio of the 
terminal velocity to the speed of light, and a determines 
how fast the velocity reaches a terminal value. Within 
the notation of Q , po is a function of the luminosity of 

the galaxy given by po = — 3 — , where 

7(L) 2 = 0.72 + 0.44 log L/L*. In our present work we will 
assume that both /3q and a are simply constants that can 



(x 2 + /3$x 2 + a 2 )iW 



{x 2 + 2j3 2 x 2 + a 2 ) (x 2 + a 2 ) (1 " /3o_) 

(15) 

which can be inserted immediately in Eq. and the 
equation of state can be computed. Indeed, once p(x) is 
computed, the density will be in this case given by 
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which is a decreasing function. If we demand that dark 
matter is regular matter that fulfills the positive energy 
condition, then it is sufficient to demand p(0) > 0. Hence 
Eq. (|16|) implies an upper limit on the central pressure 
completely determined by the terminal velocity po and 
a, which explicitly is: 



P(0) < 



(17) 



In order to advance with our method, we will use obser- 
vational data for the rotational curve profile observed in 
several low surface brightness galaxies, reported in (l3| . 
where for computing the velocity due to the dark halo, 
Whalo, we are considering w ha io 2 = Vobservcd 2 -Wdisk 2 , where 
^observed is the magnitude of the observed velocity, and 
fdisk is the magnitude of the velocity computed due to 
the observed luminous distribution. 

The procedure for a given galaxy, consists in determin- 
ing the best value of the parameters po and a that fits 
the observed rotational curve. Those values are reported 
in table U 

Once (3 2 (x) has been fixed for each galaxy, then solve 
the integral Eq. (|SJ). 

It is worth to notice that, for the (usual) case when 
po <C 1, the Eq. © can be readily solved: 
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We use this approximate expression to evaluate the 
initial pressure p(0) in the full description, such as p(0) = 
Papprox(O) = /3p/3ga 2 . Hence the integration constant 
can be directly computed as 



C = p(0)a 



3q 



(19) 





Galaxy 


Pa KT 4 


a 


p(0)10- 4 


p(0) 10" 12 


U D M 10" 9 


1 


UGC3137 


3.68 


4.09 


162.28 


733.98 


45.42 


2 


UGC711 


3.33 


6.62 


50.72 


188.00 


36.58 


3 


DD0189 


1.95 


2.71 


103.08 


130.69 


12.86 


4 


NGC100 


3.11 


2.68 


269.18 


868.53 


32.73 


5 


NGC2366 


1.22 


1.17 


218.71 


109.15 


5.11 


6 


NGC3247 


2.25 


1.48 


469.59 


812.25 


17.68 


7 


NGC5023 


2.97 


1.63 


662.89 


1953.05 


30.09 


8 


NGC4455 


2.58 


6.21 


34.57 


76.73 


21.96 


9 


UGC10310 


2.82 


4.03 


97.57 


258.25 


26.59 


10 


UGC1230 


2.83 


3.89 


106.02 


283.71 


26.91 


11 


UGC1281 


2.33 


3.17 


108.54 


197.01 


18.35 


12 


UGC3371 


2.85 


8.14 


24.53 


66.45 


26.49 


13 


UGC4173 


1.84 


12.46 


4.39 


5.01 


10.88 


14 


UGC4325 


5.11 


6.27 


133.26 


1163.15 


86.35 


15 


UGC5750 


3.16 


13.49 


10.98 


36.61 


31.73 


16 


LSBF563 


2.88 


4.22 


93.29 


258.56 


27.80 



TABLE I: Fits for several galaxies, with q = 510 -7 , which 
corresponds to a system with distance in Kiloparsec units and 
with a halo mass at ten Kiloparsecs of 10 10 Solar masses. 



Our results for the mass density, the pressure and the 
resulting pvsp relation are shown in figure [1] for all LSB 
galaxies listed in Table |TJ In the bottom right panel of Fig. 
[1] we have included in shadow grey region the cosmolog- 
ical bounds obtained by Muller in [3j for a barotropic 
equation of state. It is interesting that the dark matter 
equation of state is fully determined by f3 2 {x) and the 
election of p(0). If we use p(0) given by Eq. (fl"9)) and 
as the density profile Eq. (|T4"|) . the resulting equation of 
state for dark matter are compatible with the cosmologi- 
cal bounds as can bee seen in the figure. Furthermore, to 
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FIG. 1: Top left panel: Velocity profiles for each of the 16 
LSB galaxies listed in table [I] Top right panel: Pressure pro- 
files. Bottom left panel: Density profiles. Bottom right panel: 
Pressure as a function of the density, gray area correspond to 
cosmological bounds on P(p) [3]. 



first order there is a p = lojjmP functional dependence of 
the pressure as a function of the density as can be seen in 
the same figure. The resulting values of ojdm fitted from 
our the exact solutions for p are reported in last column 
of Table |U 

We can ask ourselves how large can be the pressure at 
the center of the galaxy such as 

1. p(x) > 0, i.e. fulfills positive energy condition and 

2. -1.5 x 10- e p{x) < p(x) < 1.13 x 10- e p{x) i.e. ful- 
fills cosmological bounds [|| . 

In Fig. [2] it is shown the pressure at the center of galaxy 
such that both conditions are considered. 

Finally, once the solutions for the equation of state 
within general relativity were obtained, we can study the 
Newtonian limit in order to asses its accuracy. In this 
case, given a velocity profile (3 2 (x) Eq. (jT4j) , one obtains 
directly the density and the mass function, and it is left, 
as before, to perform an integration for obtaining the 
pressure. 

It is obtained the following expressions for the mass 
function, density and pressure analytically from Eqs. 

mm- 



(20) 
(21) 
(22) 

In the Newtonian case, it is possible to obtain an analyt- 
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FIG. 2: Values for the central pressure for some typical dark 
matter dominated galaxies. Circles correspond to the election 
for the central pressure given by p(0) — /?o/3ga 2 . The red line 
shows the allowed region for the central pressure such as the 
resulting EOS satisfy both the positive energy condition and 
the cosmological bounds on barotropic equation of state found 
in [|. 



ical expression for the equation of state for dark matter 
P 2 p (fi 2 + 6a 2 qp + /? ov //?o+24aV) 



A, + v^o + 24 «V 



(23) 



In excellent agreement with the general relativistic re- 
sults. 



IV. CONCLUSIONS 

We have modeled dark matter as a perfect fluid with- 
out assuming any specific form of the equation of state. 
We have shown that, in a spherically symmetric space, 
the rotation profile curve directly implies the pressure as 
a function of the density. By demanding dark matter to 
fulfill positive energy condition, a maximum pressure at 
the center of galaxies was obtained, namely: 



p(0) < 1.22 x 10 



-16 



1 grms 



(24) 



where Vq is the asymptotic circular velocity of the halo. 

Furthermore, we have shown that for certain values of 
p(0), the resulting equation of state is consistent with 
cosmological bounds for a barotropic equation of state 
(see Fig. [2]) . Departure from the pressure-less hypothesis 
(p 0) for the dark matter at the astrophysical and 
cosmological level may help us to understand the nature 
of the dark matter. 

In this work, we shown that the observations make 
plausible that the dark fluid has a relation between the 
pressure and the density of the form p = ujjjm p, where 
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udai was obtained from a linear fit to the EOS obtained 
from each galaxy (See last column of Table|T]for numerical 
values). Actually, the density, p must be expressed in 
terms of the rest mass density, po, the one which satisfies 
the continuity equation, {pqu 11 )-^ = 0, and the internal 
energy of the fluid, e, obtaining: 



V = Udm Pqc 2 (1 + e). 



(25) 



If we deal with the dark fluid as if it were a perfect fluid, 
then we obtain an expression which allows us to relate the 
velocity parameter, /3q 2 (identified with ujdm), in terms 
of well defined thermodynamical functions: 



RT 



mo c 2 (1 + e) 



A) 



(26) 



so that we can understand the dark matter halo as formed 
of a dark fluid with a temperature. The parameters of 



the observed rotational velocity profile are related to such 
temperature, so that there is not a different type of dark 
fluid for each halo, but the same dark fluid at different 
temperatures. 

We will give more ground to this lucubration in a forth- 
coming work. 
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